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CRESTED PRODUCTS OF MARKOV CHAINS 

By Daniele D'Angeli and Alfredo Donno 

Universite de Geneve 

In this work we define two kinds of crested product for reversible 
Markov chains, which naturally appear as a generalization of the case 
of crossed and nested product, as in association schemes theory, even 
if we do a construction that seems to be more general and simple. 
Although the crossed and nested product are inspired by the study 
of Gelfand pairs associated with the direct and the wreath product 
of two groups, the crested products are a more general construction, 
independent from the Gelfand pairs theory, for which a complete 
spectral theory is developed. Moreover, the fc-step transition prob- 
ability is given. It is remarkable that these Markov chains describe 
some classical models (Ehrenfest diffusion model, Bernoulli-Laplace 
diffusion model, exclusion model) and give some generalization of 
them. 

As a particular case of nested product, one gets the classical Insect 
Markov chain on the ultrametric space. Finally, in the context of the 
second crested product, we present a generalization of this Markov 
chain to the case of many insects and give the corresponding spectral 
decomposition. 

1. Introduction. The starting point for this work is the paper [4] on 
Gelfand pairs (for a general theory and applications see [7]) by Ceccherini- 
Silberstein, Scarabotti and Tolli. In that paper, particular cases of homo- 
geneous spaces are presented. Given two finite groups G and F, together 
with two subgroups K < G and H < F, let X = G/K and Y = F/H be the 
corresponding homogeneous spaces. The Gelfand pairs associated with the 
direct product G x F and the wreath product F I G are studied. A complete 
spectral decomposition, together with the associated spherical functions, is 
given for the space L{X x Y). 

One can ask if it is possible to get the same spectral decomposition with- 
out using the group structure, but considering the product (in some sense) 
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of two reversible Markov chains P and Q defined on two finite sets X and 
y, respectively. 

The analogue of direct product and wreath product is developed in Sec- 
tions 5 and 6, respectively. In particular, in Section 6, the Markov chain on 
the ultrametric space, studied by Figa-Talamanca in [II] and that we call 
"Insect Markov chain," is revisited from a different point of view. See also [6]. 

An equivalent approach to Gelfand pairs product comes from the asso- 
ciation schemes theory (see [2]). In particular, the direct and the wreath 
product of groups correspond to the crossed and nested product of associa- 
tion schemes. In [3] Bailey and Cameron introduced a more general product 
for association schemes called the crested product. 

Inspired by this definition, in this paper two new products of Markov 
chains are introduced. We will refer to them as the first and the second 
crested product, respectively. Actually, the definition does not seem to be 
the exact analogue, in terms of Markov chains, of the crested product for 
association schemes, but it looks even more interesting. In fact, they ex- 
tend the crossed and the nested product and a complete spectral theory is 
presented in Sections 4 and 7. 

In particular, in Section 7, we show that the second crested product is 
a generalization of a Markov chain that we will call multi-insect, which is 
an analogue of the so-called nonbinary Johnson scheme, already studied in 
[4], that corresponds to considering h insects living in different subtrees 
and moving only one per each step in such a way that their distance is 
preserved. We give a generalization of the results of [4] in Section 7 with 
a more elementary and direct proof. This problem has been suggested by 
Scarabotti. As an example, the case of the bi-insect is studied. 

2. Basic examples and motivations. In this section we present some ex- 
amples, classical and new, of finite diffusion processes that one can treat 
using methods from harmonic analysis. In the next sections, we will study 
them by introducing special products of finite Markov chains and performing 
a complete spectral analysis of the associated Markov operator. 

2.1. Classical models. 

2.1.1. The Ehrenfest diffusion model. The classical Ehrenfest diffusion 
model consists of two urns numbered 0, 1 and n balls numbered 1, . . . ,n. A 
configuration is given by a placement of the balls into the urns. Therefore 
there are 2" configurations and each of them can be identified with the subset 
A of {1, . . . , n} corresponding to the set of balls contained, for instance, in the 
urn 0. Note that there is no ordering inside the urns. The initial configuration 
corresponds to the situation when all balls are inside the urn 0. 
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Then, at each step, a ball is randomly chosen (with probability 1/n) and it 
is moved to the other urn. Denoting by Qn the set of all subsets of {1, . . . , n}, 
we can describe this random process as follows. If we are in a configuration 
A € Qn, then, at the next step, we are in a new configuration B € Qn of the 
form ^U{i} foi' some j ^ A or A \ {j'} for some j' G A and each of these 
configurations is chosen with probability 1/n. 

In [5], the authors give a generalization of this model to the case of m 
urns and n balls. 

2.1.2. The Bernoulli-Laplace diffusion model. The Bernoulli-Laplace dif- 
fusion model consists of two urns numbered 0, 1 and 2n balls, numbered 
l,...,2n. A configuration is a placement of the balls into the two urns, n 
balls each. Therefore there are (^") configurations; each of them can be 
identified with an n-subset A of {1, ... , 2n} corresponding to the set of balls 
contained in the urn 0. The initial configuration corresponds to the situation 
when the balls contained in the urn are 1,2, ... ,n. Denoting by tbe 
set of all n-subsets of {1, 2, . . . , 2n}, we can describe this random process as 
follows. If we are in a configuration A G then at the next step we are in 
a new configuration B G 0„ of the form j4JJ{?} \ {j}, for some i ^ A and 
j A and each of these configurations is chosen with probability 1/n^ . 

2.1.3. The exclusion model. Let {X, E) be an undirected connected graph 
with n vertices. Fix h <n and choose an /i-subset A of X. We put a particle 
on each vertex of A. Then the exclusion process is defined as a Markov chain 
on the set of /i-subsets of X. Starting from the state A, we pick a par- 
ticle on some vertex of A with probability proportional to the degree of the 
vertex that it occupies, pick a neighboring site of this vertex at random and 
move the particle to the neighboring site provided this site is unoccupied. 
If the site is occupied by another particle, the chain stays at A. See [8] and 
[12] for more examples and details. 

2.1.4. The Insect Markov chain. This Markov chain, in its classical ver- 
sion introduced by Figa-Talamanca [11], is performed on the (finite) ultra- 
metric space given by the nth level of a g-ary rooted tree. At the starting 
instant the insect stays in a fixed vertex of the nth level. It performs a simple 
random walk on the tree. The next state of the Markov chain is given when 
the insect reaches again the nth level. The probability transition between 
two vertices only depends on their ultrametric distance and so this Markov 
chain is invariant under the action of the full automorphisms group of the 
rooted tree. 
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2.2. Generalizations. 

2.2.1. The first crested product. The Ehrenfest model, generalized to the 
case of n balls and m urns, is described by the crossed product of n Markov 
chains defined on a space of cardinality m (the space of the urns), as will 
be shown in Section 5. Actually, the definition of crossed product that we 
give admits that the choice of the ith. ball is weighted by a not necessarily 
homogeneous probability distribution p^. 

One can ask what happens if a hierarchy in the set of balls is introduced. 
So suppose to have n ordered balls, numbered by 1, . . . ,n, and m urns. 

At each step, one ball (say the ith one) is randomly chosen and it is moved 
into another urn following a transition probability Pj defined on the space 
of the urns. For all j > i, the ball j is also randomly moved to another urn. 

This model is described by the nested product of n Markov chains defined 
on a space of cardinality m (the space of the urns) , as will be shown in Section 
6, where we suppose that each ball is chosen according with a probability 
distribution p^. 

The Insect Markov chain described in Section 6.3 is obtained from this 
model for a particular value of p^ and forcing the transition probability Pi 
to be uniform, as shown in Proposition 6.4. 

The first crested product introduced in Section 4 is a generalization of 
both crossed and nested product and it describes the following diffusion 
model. 

The (C, N) -Ehrenfest model. We have n balls numbered by 1, . . . , n and m 
urns. Suppose that we have a partition of the set of the balls in two subsets 
C]JN. The balls in C are white, the balls in are black. At each step, we 
choose a ball i according with a probability distribution p^: if it is white, 
then we move it to another urn following a transition probability Pi and all 
the other balls are not moved. If it is black, we move it by Pi and then each 
ball (white or black) numbered by j, for j > i, is moved uniformly to a new 
urn. 

The card players model. The previous model can also be described in the 
following way. Suppose that n card players sit at the same edge of a table. 
Each of them has a deck of cards consisting of m cards. Suppose that the 
set of players is partitioned in two subsets: blond hair and brown hair. At 
the starting moment, each player has a card. A player is randomly chosen, 
so he puts his card inside the deck and then chooses another card. If he is 
blond, the others keep their own card; if he is brown, then all the players 
(blond or brown) sitting at his right side also must change their card. 
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2.2.2. The second crested product. 

The bi-insect. This is a generalization of the Insect Markov chain. In this 
case we have two insects that, at the starting time, hve in the nth level of 
a tree, following the condition that their (ultrametric) distance is maximal, 
which is equivalent to requiring that they live in two different subtrees of 
depth n — 1. At each step, with probability pq they do not leave their own 
subtrees; with probability 1 — Po one of them can leave its subtree choosing 
another subtree not occupied by the other insect. In the first case, we ran- 
domly choose one insect and change its position according with the Insect 
Markov chain on the corresponding subtree. In the second case, the moved 
insect randomly chooses its new position in the new subtree that it has occu- 
pied. This yields a Markov chain on the space of all possible configurations 
of two insects at maximal distance in a tree of depth n; this Markov chain 
will be studied in Section 7.2. 

The multi-insect. An immediate generalization of the previous construction 
can be performed considering h insects living in the nth level of the tree and 
having maximal distance from each other. At each step, with probability pQ 
only one insect moves to another site staying in the same subtree; with 
probability 1 — po only one insect occupies a new subtree of depth n — 1 
and randomly chooses its position. This situation will be referred to as the 
multi-insect in Section 7. 

Moreover, the multi-insect can be also regarded as an exclusion process 
on the ultrametric space with the constraint that the h insects keep maximal 
distance from each other. At each step an insect moves, as in the exclusion 
process a particle does. 

The aim of this work is also to diagonalize a more general Markov chain, 
that we call the second crested product, described in Section 7 and that 
corresponds to the following model. 

The crested Bernoulli-Laplace model. Suppose that we have two urns num- 
bered 0, 1 and n balls: the urn contains h balls and the urn 1 contains the 
remaining n — h balls. The h balls in the urn are endowed with a label 
that can go from 1 to s (note that two different balls can have the same 
label). Moreover, we can use a (nonhomogeneous) coin. At each step, we flip 
the coin: if we get heads, then we randomly choose a ball i in the urn and 
a ball j in the urn 1 and exchange them. Before putting j in the urn we 
attach a label (randomly chosen in {1, . . . , s}) on it. 

If we get tails, then we randomly choose a ball i in the urn and we 
change or not its label according with a transition probability Q on the set 

Observe that, if we get heads, then we perform on the set of the balls an 
analogue of the Bernoulli-Laplace diffusion model, where the urns contain 
h and n — h balls, respectively. 
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Although in Section 7.1 the spectral analysis of the second crested product 
is given, in Section 7.2 we give explicit computations only in the case of the 
bi-insect, which seems to be easier to handle. 

2.3. Motivations. The classical models described in the previous sections 
are well known in literature and they have been studied by many authors 
using approaches deriving from Gelfand pairs theory, representations theory, 
associations schemes, distance regular graphs and orthogonal polynomials. 

In the Ehrenfest model, generalized to the case of m urns and n balls, 
the set of all configurations can be identified with the Hamming scheme 
Xn,m = {0, ...,m — 1}" consisting of all the functions defined on the set 
{!,..., n} taking values in {0,...,m — 1}. The group that naturally acts 
on Xn,m is the wreath product Sm I Sn- The stabilizer of a fixed function 
in Xn^m is isomorphic to the subgroup 5^-1 I Sn- In [5], for example, the 
authors show that {Sm I Sn, Sm-i I Sn) is a Gelfand pair, they give the de- 
composition of the space L{Xn,m) of the complex functions defined on Xn,m 
into irreducible submodules and provide an expression for the respective 
spherical functions by using Krawtchouk polynomials. These submodules 
coincide with the eigenspaces of the Markov operator on L(X„ m) describ- 
ing the Ehrenfest model. A different approach to the Hamming scheme is 
given in [2] by using association schemes theory. 

Following the idea of Diaconis, it is possible to study the rate of conver- 
gence to the stationary distribution of finite Markov chains by using Gelfand 
pairs theory. In particular, given a Markov chain which is invariant under the 
action of a group G, its transition operator can be expressed via a Fourier 
series where the classical exponentials are replaced by the irreducible repre- 
sentations of G. In [10] the authors apply this argument to the generalized 
Bernoulli-Laplace diffusion model, consisting of two urns containing h and 
n — h balls, respectively. All the states of this Markov chain constitute the 
so-called Johnson scheme associated with the Gelfand pair {Sn,Sh x Sn-h)- 
The corresponding decomposition into irreducible subrepresentations (and 
so into eigenspaces for the associated Markov operator) is given, for instance, 
in [5], where the spherical functions are described in terms of Hahn polyno- 
mials. Like the Hamming scheme, also the Johnson scheme can be presented 
as an association scheme (see [2]). 

Finally, the Insect Markov chain can be studied using a Gelfand pair 
associated with the automorphisms group of the rooted tree (see [6] and 
[11]). The Fourier analysis that one can develop allows to give a complete 
spectral analysis and to investigate some ergodic properties. 

The approach proposed in this paper is more elementary. We define and 
study several Markov chains, generalizing the classical models described 
above, by using only techniques and arguments coming from linear and 
multilinear algebra. This allows to give directly the spectral analysis of the 
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associated Markov operator, without invoking the theory of groups and their 
representations. Our construction shows what reahy are the structures to the 
base of the analysis of these models and suggests the way to approach more 
general cases. 

In the case of the first crested product we consider the Cartesian product 
X = Xi X • • • X Xn of n finite spaces and a Markov operator P on L{X) 
given by a convex combination of Markov operators on L{Xi) depending 
on a partition of the set {1, . . . ,n}. The eigenspaces of the operator P are 
expressed as certain combinations of the tensor products of the eigenspaces 
of the Markov operators on every space L(Xj). Their form is easier in the 
case of the crossed product when none of the spaces Xi is distinguished 
from the other ones (see Section 5). On the other hand, it becomes more 
complicated if one introduces a hierarchy. The extremal case gives the nested 
product (see Section 6), where the space X can be identified with the nth 
level of a rooted tree. 

In the second crested product we consider the Cartesian product of two 
finite spaces X and Y , where we are given the eigenspaces of a Markov 
operator Q defined on L{Y). Set 0/j the space of the functions with domain 
an /i-subset of X and values in Y . We define a Markov operator P [see 
(6)] on L{Qfi) and we investigate its spectral decomposition by using two 
differential operators D and D* (see Definition 7.1) and the known spectral 
analysis of Q. 

Although the computational aspect can be more or less complicated de- 
pending on the singular cases, the methods used do not require anything 
more than linear or multilinear algebra. 

3. Preliminaries. The following topics about finite Markov chains can 
be found in [6]. 

Consider a finite set X, with \X\ = m. Let P be a stochastic matrix of 
size m whose rows and columns are indexed by the elements of X, so that 

^ p{xQ,x) = 1, 

for every xq ^ X. Consider the Markov chain on X with transition matrix 
P. 

Definition 3.1. The Markov chain P is reversible if there exists a strict 
probability measure vr on X such that 

7r(x)p(x,y) =TT{y)p{y,x), 

for all x,y G X. 



8 D. D'ANGELI AND A. DONNO 

We will say that P and vr are in detailed balance. For a complete treatment 
about these and related topics see [1]. 

Define on L{X) = {f -.X — > C} a scalar product in the following way: 

(/i, 72)77 = /i(2;)/2(a:)vr(x), 
xex 

for all /i, /2 G L{X) and the linear operator P : L{X) — > L{X) by 

{Pf){x)=Y.p{x,y)f{y). 

y&X 

It is easy to verify that vr and P are in detailed balance if and only if P is self- 
adjoint with respect to the scalar product {■,-)n- Under these hypotheses, it 
is known that the matrix P can be diagonalized over the reals. Moreover, 1 
is always an eigenvalue of P and for any eigenvalue A one has |A| < 1. 

Let \z be the eigenvalues of the matrix P, for every z X , with A^q = 1. 
Then there exists an invertible unitary real matrix U = {u{x,y))x^y£x such 
that PU = UA, where A = {XxSx{y))x,yex is the diagonal matrix whose 
entries are the eigenvalues of P. This equation gives, for all x,z ^X, 

(1) ^ ^(2;, =u(x,2;)A^. 

yex 

Moreover, we have U^DU = I, where D = {'K{x)6x{y))x,yGX is the diagonal 
matrix of coefficients of vr. This second equation gives, for all y,z (zX, 

(2) 'Yu{x,y)u{x,z)7r{x) = 6y{z). 
xex 

Hence, the first equation tells us that each column of U is an eigenvector of 
P; the second one tells us that these columns are orthogonal with respect 
to the product 

Proposition 3.2. The kth step transition probability is given by 

(3) p^''\x,y) = Tr{y) ^ u{x , z) X';u{y , z) , 

zex 

for all x,y X . 

Proof. The proof is a consequence of (1) and (2). In fact, the matrix 
U^D is the inverse of U, so that UU^D = I. In formute, we have 

V u{x,y)u{z,y) = -^A^(x). 
^x ^(^) 

From the equation PU = UA we get P = UAU'^D, which gives 

p{x,y)=TT{y) u{x,z)Xzu{y,z). 
zex 



CRESTED PRODUCTS OF MARKOV CHAINS 



9 



Iterating this argument, we get 

P'' = UA'^U^D, 

which is the assertion. □ 

Recall that there exists a correspondence between reversible Markov chains 
and weighted graphs. 

Definition 3.3. A weight on a graph Q = {X, E) is a function w.X x 
X — > [0, +oo) such that: 

(1) w{x,y) = w{y,x); 

(2) w{x, y) > if and only if x ~ y. 

If ^ is a weighted graph, it is possible to associate with w a stochastic 
matrix P = {P{x,y))x^y£x on X by setting 

. X w{x,y) 

with VF(a;) = X^^ex '"^(^' -^)- "^^^ corresponding Markov chain is called the 
random walk on Q. It is easy to prove that the matrix P is in detailed 
balance with the distribution vr defined, for every x G X, as 

Wix) 

TTix) — 



W ' 

with W = J2z£X W{z). Moreover, vr is strictly positive if X does not contain 
isolated vertices. The inverse construction can be done. So, if we have a 
transition matrix P on X which is in detailed balance with the probability 
TT, then we can define a weight w as w{x,y) = Tr{x)p{x,y). This definition 
guarantees the symmetry of w and, by setting E = {{x, y} : w{x, y) > 0}, we 
get a weighted graph. 

There are some important relations between the weighted graph associ- 
ated with a transition matrix P and its spectrum. In fact, it is easy to prove 
that the multiplicity of the eigenvalue 1 of -P equals the number of connected 
components of Q. Moreover, the following propositions hold. 

Proposition 3.4. Let Q = {X, E, w) be a finite connected weighted graph 
and denote P the corresponding transition matrix. Then the following are 
equivalent: 

(1) Q is bipartite; 

(2) the spectrum (t{P) is symmetric; 

(3) -iea(P). 
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Definition 3.5. Let P be a stochastic matrix. P is ergodic if there 
exists no € N such that 

p^"") (x, y) > for ah x,yeX. 

Proposition 3.6. Let Q = {X,E) be a finite graph. Then the following 
conditions are equivalent: 

(1) Q is connected and not bipartite; 

(2) for every weight function on X, the associated transition matrix P is 
ergodic. 

So we can conckide that a reversible transition matrix P is ergodic if and 
only if the eigenvalue 1 has multiplicity 1 and —1 is not an eigenvalue. 

In what follows we always suppose that the eigenvalue 1 has multiplicity 
1, so that the graph associated with the probability P is connected. This is 
equivalent to requiring that the probability P is irreducible, according with 
the following definition. 

Definition 3.7. A stochastic matrix P on a set X is irreducible if, for 
every xi,X2 €X, there exists n = n{xi,X2) such that p^'^\xi,X2) > 0. 

4. The first crested product. In this section we introduce a particular 
product of Markov chains defined on different sets. This idea is inspired by 
the definition of crested product for association schemes given in [3] . 

Let Xi be a finite set, with |Xj| = mj, for every z = 1, . . . , n, so that we can 
identify Xi with the set {0, 1, . . . , mj — 1}. Let Pj be an irreducible Markov 
chain on Xi and let pi be the transition probability associated with Pj. 
Moreover, assume that pi is in detailed balance with the strict probability 
measure on Xi, that is, 

ai{x)pi{x,y) =ai{y)pi{y,x), 

for all x,y € Xi. 

Consider the product XiX ■ ■ - x X^. Let {1, . . . ,n} = C]JN he a partition 
of the set {l,...,re} and let Pi,P2: ■ ■ ■ ,Pn be a probability distribution on 
{1, . . . ,n}, that is, Pi>0 for every i = 1, . . . ,n and J27=iPi = 1- 

Definition 4.1. The first crested product of Markov chains Pj's with 
respect to the partition {1, . . . ,n} = CIJN is the Markov chain on the prod- 
uct Xi X • • • X Xn whose transition matrix is 

P = "^P^ih • • • ® /j-i Pj li+i (8) • • • (8) /„) 

iec 

+ p^iih ® • • • ® (g) Pi ® Ji+i (g) • • • J„), 
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where Ij denotes the identity matrix of size nii and Ji denotes the uniform 
matrix on Xi, that is, 



Ji 



1 

mi 



1 ■•. 

Vi ... 



1\ 



1/ 



In other words, we choose an index i in {1, . . . ,n} following the distribu- 
tion . . . If i G C, then P acts on the ith coordinate by the matrix 
Pi and fixes the remaining coordinates; if i G A^, then P fixes the coordi- 
nates — 1}, acts on the iih. coordinate by the matrix Pi and changes 
uniformly the remaining ones. 

For all (xi, . . . , Xn), {yi, ■ ■ ■ , Un) G Xi x ■ ■ ■ x Xn, the transition probability 
p associated with P is given by 

p{{xi,...,Xn),{yi,---,yn)) 

= ^Pi(.^i{xi,yi) ■ ■■6i-i{xi-i,yi-i) 



X Pi{xi,yi)Si+i{Xi+i,yi+i) ■ ■■Sn{xn,yn)) 
Qf5i{xi,yi) ■ ■ ■6i-i{xi-i,yi-i)pi{xi,yi) 



n 



where 5i is defined by 



Siixi,yi) 



j=i+i "^j 



if Xi = yi, 
otherwise. 



We want to study the spectral theory of the operator P. We recall that 
the following isomorphism holds: 



L(Xi X . . . X X„) 



)L{Xi), 



i=l 



with (/l (g) • • • (g) fn){xi, ...,Xn):= fl{xi)f2{x2) ' ■■fn{Xn)- 

Assume that, for every i = 1, . . . ,re, the following spectral decomposition 
holds: 

that is, Vj^ is an eigenspace for Pi with associated eigenvalue Xj^ and whose 
dimension is m,- . 

Now set N = {ii, . . . and C = {ci, . . . ,c/i}, with h + I = n and such 
that ii < ■ ■ ■ <ii and ci < ■ ■ ■ < Ch- 
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Theorem 4.2. The probability P defined above is reversible if and only 
if Pk is symmetric for every k> ii. If this is the case, P is in detailed balance 
with the strict probability measure vr on Xi x • • • x Xn given by 

T^\X\ , . . . , Xn ) — 

rriii+i ■■■mn 

Proof. Consider the elements x = (xi, . . . , Xn) and y = {yi, . . . , yn) be- 
longing to Xi X • • • X Xn- First, we want to prove that the condition cjfc = — , 
for every k > ii, is sufficient. Let A: € {1, . . . , n} such that Xi = yi for every 
i = 1, . . . , fc — 1 and Xk 7^ yk- Suppose k <ii. Then we have 

p{x, y) = Pk{Pk{xk,yk)Sk+i{xk+i,yk+i) ■ ■ ■ K{Xn, Vn))- 
If Xi = yi for every z = /c + 1, . . . , n, we get 

iT{x)p{x,y) = ai{xi) ■ ■■ak{xk) ■ • • o-j^(xi Jp^.- 



Pk{xk-,yk) 



t ^ , ^ / \ Pk{yk,Xk) 
cri{yi) ■ ■ ■ ak{yk) ■ ■ ■ CTniViJPk- 



mi+i ■■■mn 
= -n{y)p{y,x), 

since ak{xk)pk{xk, Vk) = (^kiyk)Pk{yk,Xk)- If the condition Xi = yi is not satis- 
fied for every i = k + 1, . . . ,n, then the equality ■k{x)p{x, y) = iT{y)p{y, x) = 
easily follows. 

li k = ii, then we get 

p{x,y)=Pi (pi^{xi^,yi^) 



and so 



/ \ / \ / \ / \ Ph (Xii , yii ) 

7r{x)p{x, y)=ai[xi)--- ai^ {xi^ )pi — ^ 



■■■^^ 



/ \ r \ PiiKyii^Xii) 

Myi) f^ii (yii )Ph — 2 



= 'K{y)p{y,x), 

since Ui^ (x^, )pi^ {xi^,yi^) = ai^ {yi^ )pi^ {yi^,Xi^). 
In the case k > ii, we have 

p{x,y)= 2^ Pi 
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and so 



TT{x)p{x,y) = 2^ 

m ■ , -1 . . . m * ^ 



rrir. 



p 



= Tr{y)p{y,x). 

In fact, the terms corresponding to an index i < k satisfy pi{xi,yi) = Pi{yi,Xi) 
since Xi = yi; the term corresponding to the index k satisfies Pk{xk-,yk) = 
PkiUkiXk) since the equahty 

Pk{xk,yk) =Pk{yk,Xk) 

holds by hypothesis. 

Now we want to prove that the condition ak = — , for every k > ii, is 
also necessary. Suppose that the equality ■K{x)p{x,y) = 7r(y)p(y, x) holds. 
By the hypothesis of irreducibility we can consider two elements x^,y^ G 
XiX ■ ■ ■ X Xn such that x^^ ^ y^^ and with the property that pi^ (x?^ iVii) 0- 
Now we have 



This gives 



7r{x')p{x',y')=7r{y^)p{y^,x') ^ 
^(^°K«,y°)=vr(2/°K(y°,<). 

Ax') _PHiyl,xl) _ai,{xl) 



Consider now the element x = {x\, . . . , x,? , y^-^j^i, ■ • ■ , y^)- The equality 7r(x)p(x, 
y^) = ir(y^)p{y^ ,x) implies 

TTjx) _ pi,{yl,x^ij _ ai.jx^ij 

niy^) ^niyV 
So we get 7r{x^) = '7r(x), that is, the probability vr does not depend on the 
coordinates + 1, . . . ,n. Set now x' = (x?, . . . ,x^^, . . . ,x^_i,Xk, ■ ■ ■ ,Xn)- The 
equality 7r(x'^)j>(x'^, x') = 7r(x')p(x', x*^) gives 

vr(x°)f E P°(p,(^°,x;.)))=vr(x')f E ^°))y 
^j(^N,j<k ) \j&N,j<k / 

Since the probability vr does not depend on the coordinates ii + 1, . . . , n, we 
get Pk{xl,x'jJ =pk{x'j^,xl). This implies cTfc(x'^) = crA;(x^) and so the hypoth- 
esis of irreducibility guarantees that ak is uniform on Xk- This completes 
the proof. □ 
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Theorem 4.3. The eigenspaces of the operator P are given by: 

W^^---(S) W^-^ Vjl V^+^ V^+^ ^■■■^V^, with jk + 0, for k G 
{ii + 1, . . . , n} and where 



with eigenvalue 



\vl, ji = 0,...,ri, ifieC, 



i€C:i<k i>k 



• Vj{ (g) ■ ■ ■ (S) vjl _l (g) vjl (g> (g) ■ ■ ■ (S) Vf^ , with jt = 0,...,rt, for every 

t=l,...,ii, with eigenvalue 

ii n 
i=l i=ii+l 

Proof. Fix an index A; G {ii + 1, ii + 2, . . . , n} and consider the function 
(p in the space 

W^®---® W^-^ ® vjl (g) Vo^+^ (8) V^+^ (8) • • • ® Fo", 
with jk 7^ and 

jL{Xi), ifieN, 
\vj^, j, = 0,...,ri, ifiGC, 



so that 93 = (g) • • • (fk-i ® ^k ^k+i (g ■ ■ ■ ® with Lpi G for i - 
1, . . . , fc — 1, G V^^ and y?; G Vq for I = k + 1, . . . , n. Set x = (xi, . . . , 
and y = {yi,...,yn); then 

(P(^)(x) = Ep(2;>y)95(y) 

3/ 

y \i€C 

X Pi(xi,yi)(5i+l(xi+i,yi+i) • ■■5nixn,yn) 

1 1 A 



+ ^Pi^i{xi,yi) ■ ■ ■ 5i-i{xi-i,yi-i)pi{xi,yi)- 



X (fiiyi) ■ ■ ■ ipk-i{yk-i)vk{yk)vk+i{yk+i) ■ ■■^n{yn) 

E \}lP^iPi{Xi^yi)'fi{yi)] 
i£C,i<k \ yi / 
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X (pi{xi) ■ • •V3i_i(xi_i)(^i+i(xi+i) • ■■ifniXn) 



ieC,i>k \ Vi ) 



+ I v^Vi^'^i.yi)— '^—^i[yi)---^n{yn 

X 991 (Xl) • ■ ■ ipi-l{Xi-l) 

1 1 

+ XN{k) PkPk{xk,yk)-^ 



X (pi{xi) ■ ■ ■ <fk~l{Xk~l)<fk{yk) ■ ■ ■ fn{yn) 
ieC,i<k i€C,i>k 

+ E [YPiPi^^i^yi^'fii^y'n 

iGN,i>k \ yi / 

X ipi{xi) ■ ■ ■ Lpi^i{xi_i)ipi+i{xi+i) ■ • • (/7„,(x„) 

+ XN{k)J2PkPk{Xk,yk) 

yk 

X (pi{xi) ■ ■ ■ (fk-l{xk-l)'Pk{yk)<fk+l{xk+l) ■ ■ ■ <fniXn) 

: Y P'^M^)+ E P'^i^) 
i&C,i<k i€C,i>k 

+ E PiV^i^)+XN{k)pl\jMx) 

ieN,i>k 



E Pi^h+Pk^jk+T.p^'pi^)^ 

Cd<k i>k / 



\i(^C,i<k i>k 

where xn is the characteristic function of N. Note that in this case the 
addends corresponding to the indices i < k, i N, are equal to since we 
have supposed jk / 0. 

Consider now the function ip in the space 

with jt = 0, . . . ,rt, for every t = 1, . . . , ii . In this case we have 

{P(p){x) =Ypix,y)^{y) 
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i£C,i<ii \ Vi / 



X (pi{xi) ■ ■ ■ (fi-i{xi_i)(pi+i{xi+i) ■ ■■(fn{Xn) 
X ipi{xi) ■ • •(/7i_i(xi_i)v9i+i(xi+i) • ■■ipn{Xn) 



ieC,i>ii \ Vi 



+ E ( E P^Pii^i^yi)— —^i{yi)---<fn{Xn)] 

X fl{xi) ■ ■ ■ ipi-l{Xi_i) 
+ E {phPhi^il^yh) ^—^iAyii)---'^n{Xn)] 



X ifiixi) ■ ■ ■ ipi^-l{Xi^-l) 

= E p'^^M^)+ E 

i(zC,i<i\ i&C,i>ii 

+ E p'ifix)+Ph^h,'Pi^) 

ieN,i>ii 

(h n \ 

T.p'K+ E p1)^^^)- 
i=l i=h+l / 

Observe that, by computing the sum of the dimensions of these eigenspaces, 
we get 

n 

mi • ■ •mfc_i(mfc - 1) +mim2 • • -mji = mim2 • ■■rrin, 

fc=ii+l 

which is just the dimension of the space Xi x ■ • ■ x Xn- □ 

Remark 4.4. The expression of the eigenvalues of P given in the pre- 
vious theorem tehs us that if Pi is ergodic for every i = 1, . . . ,n, then also 
P is ergodic, since the eigenvalue 1 is obtained with multiplicity 1 and the 
eigenvalue —1 can never be obtained. 

We can now give the matrices U^D and A associated with P. For every 
i, let C/j, Di and Aj be the matrices of eigenvectors, of the coefficients of 
CTj and of eigenvalues for the probability Pi, respectively. The expression 
of the matrix [/, whose columns are an orthonormal basis of eigenvectors 
for P, easily follows from Theorem 4.3. In order to get the diagonal matrix 
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D, whose entries are the coefficients of vr, it suffices to consider the tensor 
product of the corresponding matrices associated with the probabihty Pi, for 
every i = 1, . . . , n, as it follows from Theorem 4.2. Finally, to get the matrix 
A of eigenvalues of P it suffices to replace, in the expression of the matrix P, 
the matrix Pi by Aj and the matrix Jj by the corresponding diagonal matrix 
jf^'^^, which has the eigenvalue 1 with multiplicity 1 and the eigenvalue 
with multiplicity rm — 1. So we have the following proposition. 

Proposition 4.5. Let P be the crested product of the Markov chains 
Pi , with i = 1, . . . , n. Then we have: 



U = Efc=ii+i Ml • • • Mfc_i (Uk - Ak) (E) Ak+i <i 
• • • i8) C/ji (8) Ai+i (8) • • • (8) An, with 

( r^^»-norm if i G N 

\Ui, ifi^C, 

where 



An + Ui® U2 



tct, -norm 



1 



v 



\/cri{mi 

By Ai we denote the matrix of size rui whose entries on the first column 
are all 1 and the remaining ones are 0: 

I) = (8)^1 A. 



A,- (8 J' 



1 

diag 
i+l 



Observe that another matrix U' of eigenvectors for P is given by U' = 
^i'=i Ui. The matrix U that we have given above seems to be more useful 
whenever one wants to compute the kth step transition probability p^'^^ (0, x) 
using the formula (3), since it contains a greater number of zeros in the first 
row with respect to U' and so a small number of terms in the sum are 
nonzero. 

Suppose X = (xi, . . . , Xn) and y = (yi, . . . , y^) elements in X = Xi x ■ ■ ■ x 
Xn- From (3) and Proposition 4.5, we have 

p^^\x,y) 



-n{y) 



z&X \r=ii+l 
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X aj._|_i(Xr-(_i, S^r+l) ■ ■ ■ Oj-aiXn^ ^n) 
+ Ui{xi,Zi)- ■ ■Ui^{Xi^,Zi^) 

X aii+i(a;ii+i,2;ji+i) • ■ ■ an{xn-,Zn)^\^z 

x( ^ mi{yi,Zi)---mr-l{yr-l,Zr-l){Ur-ar){yr,Zr) 
\r=ii+l 

X ar+l{yr+l,Zr+l) ' ■■aniVnjZn) 

+ ui{yi,zi)---Ui-, {yi^.zi^) 

X aii+i(yn+i,2n+i) ' • • «n(yn, ^n) j , 

where mi,Ui,ai are the probabUities associated with the matrices Mi,Ui,Ai 
occurring in Proposition 4.5. 

5. The crossed product. The crossed product of the Markov chains Pi's 
can be obtained as a particular case of the crested product, by setting C = 
{l,...,n} (so that = 0) and it is also called the direct product. The 
analogous case for product of groups has been studied in [9] . 

In this case, we get the following transition probability: 

n 

p{{xi, . . .,Xn), (yi, . . .,yn)) = ^P°(5(xi, yi) • • ■pi{Xi,yi) ■ ■ ■d{Xn,yn)- 

1=1 

This corresponds to choosing the ith coordinate with probability and to 
changing it according with the probability transition Pj. So we get 

p((x„ . . . (y„ . . . ,y„,)) = |^°^^(-^'^^)' '[-^ ^ + ^' 

L U, otherwise. 

So, for Xi = • • • = Xn =: X and p\ = ■ ■ ■ = p%, = \ ^ the probability p defines 
an analogue of the Ehrenfest model, where n is the number of balls and 
|X| = m is the number of urns. In order to obtain a new configuration, we 
choose a ball with probability 1/n (let it be the zth ball in the urn Xj) and 
with probability Pi{xi,yi) we put it in the urn y^. 

As a consequence of Theorem 4.2, we get that if Pi is in detailed balance 
with TTj, then P is in detailed balance with the strict probability measure vr 
on Xi X • • • X Xn defined as 

7r(xi, . . . , x„) = 7ri(xi)7r2(x2) • • •7r„(x„). 

The matrix P associated with the probability p is given by 

n 

(4) P = ^p°(/i® •••(8)Pi®---(8)I„). 

i=l 
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The following proposition studies the spectral theory of the operator P 
and it is a straightforward consequence of Theorem 4.3. 

Proposition 5.1. Let (^q = 1^^, V'l; • ■ • ) V^m—i eigenf unctions of 
Pi associated with the eigenvalues Xq = 1, X\, . . . , Xl^._i, respectively. Then 
the eigenvalues of P are the mim2 ■ ■ ■ rrin numbers 

n 
k=l 

with / = (zi, . . . , i„) G {0, . . . , mi — 1} x • • • x {0, . . . , m„ — 1}. The correspond- 
ing eigenf unctions are defined as 

(fliixi, . . .,Xn)) = f},{xi) ■ ■ ■ fi^iXn)- 

As a consequence of Proposition 4.5, in order to get the matrices U,D 
and A associated with P, it suffices to consider the tensor product of the 
corresponding matrices associated with the probability Pi, for every i = 
l,...,n. For every i, let Ui, Di and Aj be the matrices of eigenvectors, of 
the coefficients of tTj and of eigenvalues for the probability Pj, respectively. 
We have the following corollary. 

Corollary 5.2. Let P he the probability defined in (4). Then we have 

PU = UA, 
U^DU = I, 

where U = ^^^^ A = {g)^^i A, and D = (g)^^i A- 

In particular, we can express the A;th step transition probability matrix 

as 

/n \/" \k/n / ^ \ 

Let x = {xi,...,Xn) and y= (yi,...,yn)- Then we get 

p^''\x, y) = TT{y) J2 ^i{x)\^ivi{y) 
I 

= 7ri(yi)---7r„(y„) 

xY^^U^i) ■ --vlMiplXi + ■ ■ ■+pl\lf^\M) ■ ■ ■ ^Uyn), 
I 

with I = {ii, . . . ,in). 

As we said in Remark 4.4, if the matrix Pi is ergodic for every i = 1, . . . , n, 
then also the matrix P is ergodic, since the eigenvalue 1 can be obtained 
only by choosing / = 0" and the eigenvalue —1 can never be obtained. 
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Remark 5.3. Put n = 1 and set X = {0,1, . . . ,m — 1}. Consider the 
action of the symmetric group Sm on X. The stabihzer of a fixed element 
xo = is isomorphic to the symmetric group Sm-i. It is well known (see [5]) 
that {Sm, Sm-i) is a Gelfand pair and the following decomposition of L{X) 
into irreducible representations holds: 

L{x) = VoeVi, 

where Vq = C is the space of constant functions on X and Vi = {f :X — > 
C : Z^i^o^ /(^) ~ 0}. So we have dim Vq = 1 and dim Vi = m — 1. 

Analogously, one can consider the action of the wreath product Sm I Sn 
on X" = X X ■ ■ ■ X X , defined in the natural way, and then one can study 
the decomposition of L{X^). We have 

n 

i=o 

with 

Vi^®Vi^®---®Vi„, 

w{ii,...,in)=j 

where w{ii, . . . ,in) = ^{k : ik = 1}. So we have dim VKj = (") (m — 1)-' . 

If we define on X the uniform transition probability, that is, Pu{x, y) = ^ 
for all x,y X, then the matrix is the matrix J of size m. 

The eigenvalues of this matrix are 1 (with multiplicity 1) and (with 
multiplicity m — 1). The corresponding eigenspaces in L{X) are, respectively, 

= C and = {/ : X ^ C : o' /« = 0}- 

This means that, by choosing the uniform transition probability on X, one 
gets again the results obtained by considering the Gelfand pair {Sm, Sm-i) ■ 

Also in the case of X" we can find again the results obtained (see [5]) 
by considering the Gelfand pair {Sm I Sn,Sm-i I Sn)- For = J we have 
Ao = 1, Ai = • ■ ■ = Xm-i = 0. Consider now the transition probability on 
X" defined in (4), with Pi = ■ ■ ■ = Pn = ^ - The eigenfunctions tpi associated 
with the eigenvalue ^ (n — j) , with j = 0, . . . ,n, are in number of (") (m — 1)-^ . 
Moreover 

j=0 

For every j = 0, . . . ,n, these functions belong to Wj and they are a basis 
for Wj. So Wj is the eigenspace associated with the eigenvalue ^{n — j). 

More generally, consider the case of any reversible transition probability p 
on X. Let Aq = 1, Ai, . . . , Afc be the distinct eigenvalues of P and Vq = C, 
Vi, . . . ,Vk the corresponding eigenspaces. We get 

L{X''')^{Vo®Vi®---(BVk)®". 



^ ^{m-iy =J2i^ ^) {m - lyi"-^ = m" = dimL(X" 
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The eigenfunctions ipi associated with the eigenvalue ^ X]f=o "^i^i^ with X]i=o '^i 
n, are 

k 
j=0 



''0 + H h rfc 



and the corresponding eigenspaces are the tensor products in which rj copies 
oiVi, for i = 0, 1, . . . , /c, appear. Moreover, the number of different eigenspaces 
is equal to the number of integer solutions of the equation 



TQ + ri^ \-rk=n, > 0, 



SO It IS ( „ ). 

The definition of multinomial coefficient as (''o+'^i^ = (''0+ guar- 

\ ro,...,rk I ro!ri!---rfe! 

antees that 

^ f ) (diml/o)''" • • • (dim Vfc)^'= = (diml/o + ■ • ■ + dim Vfc)'^ 

as we wanted. 

6. The nested product. 

6.1. General properties. The nested product of the Markov chains Pj's 
can be obtained as a particular case of the crested product, by setting N = 
{l,...,n} (so that C = 0). The term nested comes from the association 
schemes theory (see [2]). 

Consider the product 

XiX ■■■ X Xn 

and let Pi be a transition probability on Xi. We assume that pi is in detailed 
balance with the strict probability measure tTj, for all i = 1, . . . ,n. 
The formula of crested product becomes, in this case, 

n 

(5) P = "^Piih <8) • • • ® Pi <8) Ji+i (g) Ji+2 » • • • (8) Jn). 

i=l 

Theorem 4.2 tells us that P is reversible if and only if Pk is symmetric, 
for every A; > 1, that is, T^i = for every i = 2, . . . ,n. In this case, P is 
in detailed balance with the strict probability measure vr on Xi x • • • x Xn 
given by 
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So let us assume vTj to be uniform for every i = 2, . . . ,n. The transition 
probability associated with P is 

p{{xi,...,Xn),{yi,...,yn)) 

_ PiVi{xi,yi) ^ y. ■ ■ ■ ,a;j-i), {yi, . . . ,yj-i))p^jPj{xj,yj) 

77427713 ■ ■ ■ Uln ^ rUj+i ■ ■ ■ ITln 

), (yi, . . .,yn~l))PnPn{Xn,yn)- 

As we did in the case of the crossed product, we want to study the spectral 
theory of the operator P defined in (5). 
Let 

ki=0 

be the spectral decomposition of L{Xi), for all i = 1, . . . , tt, and let Aq = 1, 
A|, . . . , A* . be the distinct eigenvalues of Pi associated with these eigenspaces. 
From Theorem 4.3 we get the following proposition. 

Proposition 6.1. The eigenspaces of L{Xi X • • • X Xn) are: 

• L(Xi)(g)- • •®L(X„_i)(g) W^^, of eigenvalue p^X^^, for kn = I, ■■■ ,rn, with 
multiplicity mi • • • 7n„_i dim(VF^ ); 

. L{Xi)(S) ■ ■ ■ ®L{Xj) Wi^^^ ® ®---®W^,of eigenvalue P°+iAi+\ + 
p^_,_2 + • ■ • + Pn; = !)•••) fj+i ^"iT-d for j = 1, . . . , n — 2, with mul- 

tiplicity mi ■ ■ - rrij dim{Wl^^^ ) ; 

• Wl^0W^®---(S)W^, of eigenvalue p^XI^+pI-\ \-pl, for ki = 0,1, ... ,ri, 

with multiplicity dim(14^^^). 

Moreover, as in the general case, one can verify that, under the hypothesis 
that the operators Pi are ergodic, for i = 1, . . . , n, then also the operator P 
is ergodic. 

The application of Proposition 4.5 to the case of the nested product yields 
the following corollary. 

Corollary 6.2. Let P be the nested product of the probabilities Pi, 
with i = 1, . . . ,n. Then we have: 

• U = Ui(^A2^---^An + El=2 • ■ • {Uk - Ak) 
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6.2. k-steps transition probability. The formula that describes the tran- 
sition probabihty after k steps in the case of the nested product can be 
simphfied using the base of eigenvectors given in Corollary 6.2 and suppos- 
ing that the starting point is = (0, . . . , 0). 

From the general formula, with the usual notation, we get 



p(^)(0,2/) 



Tr{y) 



E "Y^aiiO, Zi) ■ ■ ■ 6a,_i{0, Zr-l)iUr -ar){0,Zr) 
.z£X \r=2 

X ar+l{0,Zr+l) ■ • •a„(0,z„) 

+ ui (0,^1)02(0,^2) •••a„(0,Zn) 



\r=2 



Y^aiiVl^Zi) ■ ■ ■6^^_-^{yr-l,Zr^l) {Ur - ) (?/r- , ) 
X arJ^l{yr+l-,Zr+l)---an{yn,Zn) 

+ ui(yi, ^1)02(^2, Z2)--- an{yn, Zn) 



■ Tr{y) 



.5 iSv^i(0)---v/^r-i(0) 



(Ur — ar.)(0, Zf.) 



E'^'^i(yi'0)'^'^2(y2,2;2) 



X (5cr^_i(yr-l,^r-l)('Ur " ar){yr,Zr) 
X ar+l(yr+l,^r+l) • ■■an{yn,Zn) 



+ «i(o,^i) p?A^, + E^'m Wi(yi,^i) 

21^0 \ m=2 / 

2i=0,i>l 
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Observe that in this case the sum consists of no more than 

n n 

\X^\ + Y,{\X,\-l)=Y,\X,\-n + l 

i=2 i=l 

nonzero terms. 

Example 6.3. We want to express the feth step transition probabihty 
in the case n = 2. So consider the product X xY , with X = {0, 1, . . . , m} 
and y = {0, 1, . . . , n}. Let 

r s 

L{X) = @Vj and L{Y) = @W, 

j=0 i=0 

be the spectral decomposition of the spaces L{X) and L(Y), respectively. 
Let Ao = 1, Ai, . . . , Ar and /io = 1, /Ui, . . . , ^Ug be the distinct eigenvalues of 
Px and Py, respectively. Then the eigenspaces of L{X x Y) are L{X) (g) Wi, 
for i = 1, . . . , s, with dimension (m + 1) dim Wj and associated eigenvalue 
PyfJ'i, and Vj (8) Wq, for j = 0, . . . ,r, with dimension dim V^- and associated 
eigenvalue Aj + Py ■ 

The expression of U becomes 



U = I--^™ ^^UY-Ay) + Ux®Ay. 

In particular, let {v^ ,vl, . . . ,vl.^^y^^, . . . ,vl, . . . ,v'^.^^y^-^} and {w°,wl,..., 

^dim(iyi)'---''^i'---''^dim(H/,)} be the eigenvectors of Px and Py, respec- 
tively, that is, they represent the columns of the matrices Ux and Uy- 

Then, the columns of the matrix U corresponding to the elements (i,0) G 
{0, . . . , m} X {0, . . . , n} are the eigenvectors f * (8) (1, • • • , 1) with eigenvalue 
p^Xi + Py- On the other hand, the columns corresponding to the elements 
G {0, . . . ,m} X {0, . . . ,n}, with j = I, . . . ,n, are the eigenvectors 

0,...,0,^i=,0,...,o) 

\/crx{i)^ J 

ith place 

whose eigenvalue is Pyfij ■ As a consequence, only m + 1 + n of these eigenvec- 
tors can be nonzero in the first coordinate, so the probability p^^\{0, 0), (x, y)) 
can be expressed as a sum of m + 1 + n nonzero terms; moreover, these terms 
become m + 1 if x 7^ 0. We have 

p(^-)((0,0),(x,y)) 



7ri{x,y))[Y,v\0y{x){p^x\+p'^ 



\k 



in N 

+ ^^^^T.^^'(^)^o{x)w^{y){p'yf^,)' 
V(^x{0)axix)p{ f 
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n+1 



■ r /dim{Vi) \ 

i=0 \ a=l / 

s / ^ dim{Wj) 



+ E y ,nw r ^ ^ «;KO)'5o(xH(2/) (p^-^.f 
j=i\Vo^-y(0)cjx(2;) ^^-^ y 

6.3. T/ie insect. It is clear that the product Xi x • • • x X„ can be regarded 
as the rooted tree T of depth n, such that the root has degree mi, each vertex 
of the first level has m2 children and in general each vertex of the ith level 
of the tree has mj+i children, for every i = 1, . . . ,n — 1. We denote the ith 
level of the tree by Lj. In this way, every vertex x G Lj can be regarded as 
a word x = xi ■ ■ ■ Xi, where xj G {0,1, . . . , mj — 1}. 

We want to show that the nested product of Markov chains is the gener- 
alization of the "insect problem" studied by Figa-Talamanca in [11]. 

Let us imagine that an insect lives in a leaf x € Ln and that it performs 
a simple random walk on the graph T starting from x. 

Then there exists a probability distribution on L„ such that, for every 
y G Ln, Hx{y) is the probability that y is the first point in L„ visited by the 
insect in the random walk. If we put p{x, y) = fJ-xiv), then we get a stochastic 
matrix P = (j)(x,y))x,y^L„- Since the random walk is ^ut(T)-invariant, we 
can suppose that the random walk starts at the leftmost vertex, that we will 
call xo = (0, . . . , 0). We recall that Aut(T) is the group of all automorphisms 
of T, given by the iterated wreath product Smn ^ 'S'm„_i I • • - l Smi ■ We want 
to study this Markov chain defined on L„. 

Set ^„ = and = 00 . . . [n — i times). For j > 0, let Uj be the prob- 
ability that the insect reaches ^j+i given that is reached at least once. 
This definition implies ao = 1 and ai = ^ . In fact, with probability 1, 
the insect reaches the vertex at the first step and, starting from ^i, with 
probability ^ it reaches ^2, while with probability it returns to 

Ln- Finally, we have = 0. For 1 < j < n, there is the following recursive 
relation: 

1 , mn+i-j 
aj = h aj-iaj — . 

rUn+i-j + 1 rUn+i-j + 1 

In fact, starting at ^j, with probability ^ the insect reaches in one 

step ^j+i, otherwise with probability it reaches ij-i or one of its 

brothers; then, with probability Oj-i it reaches again and one starts the 
recursive argument. 

The solution, for 1 < j < n — 1, is given by 

l + mn + ninmn-l + mnmn-imn-2 H h m„m„_im„_2 • ■ ■ mn-j+2 

aj = — 

l + mn + mnmn-i + mnmn-imn-2 H h mn'mn-imn-2 ■ ■ ■ rUn-j+l 
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mnmn-imn-2 ■ ■ ■ rUn-j+l 



1 



1 + m„ + m„,m„_i + mnmn-imn-2 H h mnmn~imn-2 ■ ■ ■ rrin-j+i 

Moreover, we have 

p{xo,xo) = —(1 - ai) H ai(l - 02) H 

m„ mnVrin-i 

H aia2 ■ • ■an-2(l - an-ij H • • •an-l• 
m„m„_l•••m2 m„---mi 

Indeed the jth summand is the probabihty of returning back to xq if the 
corresponding random walk in T reaches but not Cj+i- It is not difficult 
to compute p{xq,x), where x is a point at distance j from xq. For j = 1, we 
clearly have p{xq,xq) =p{xq,x). We observe that, for j > 1, to reach x one 
is forced to first reach ^j, so that we have 

p(xo, x) = aia2 ■ ■ ■ aj-i{\ — a,) + • ■ • 

mn • • • rUn-j+l 

H aia2 • • •a„_2(l - a„„ij H 0102 • --an-i- 

m„ • • • m2 mn---mi 

Since the random walk is invariant with respect to the action of Aut(T), 
which acts isometrically on the tree, we get the same formula for any pair 
of vertices x,y € such that d{x, y) = j. 

Proposition 6.4. The stochastic matrix 
p{{xi,...,xn),{yi,...,yn)) 

^ PiPiixi,yi) ^ y: ^((xi, . . . , Xj-i), jyi,..., yj-.i))p'^jpj{xj,yj) 

), (yi, . . .,yn-l))PnPn{Xn,yn), 

defined in (5), gives rise to the Insect Markov chain on Ln, regarded as 
Xi X ■ ■ ■ X Xn, choosing p^ = 0102 • • •an-j(l — a„_j+i) for i = 1, . . . ,?i — 1 
and p^ = l — ai and the transitions probabilities pj 's to be uniform for all 
j = l,...,n. 

Proof. Set, for every i = 1, . . . ,n — 1, 

p° = aia2 ■ ■ ■ an-j(l - "n-i+i) 

and p^ = l — ai- Moreover, assume that the probability pi on Xj is uniform, 
that is, 



Pi — Ji' 
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If d{xo,x) = n, then we get 

Qia2---a„_i 
p{xo,x) = . 

777-17712 • • • 777„ 

If d{xo,x) = j > 1, that is, x^ = Xi for all i = 1, . . . , 77 — j, then 

p{xQ,x) = h 2^ . 

777 1 7772 •• -W-n ^ 777,^ • ■ • mj+2"li+l 

Finally, \i x = xq^ we get 

, , aia2---an~i , ^ ai • • • a„_j_i(l - a^-i) , (1-ai) 
p(xo,xo) = h > \ . 

17111712- ■■mn ^ r77„ • • • 777i+2"ii+l "7„ 

This completes the proof. □ 

The decomposition of the space L[Ln) = L[Xi x • • • x Xn) under the 
action of Aut{T) is known (see [6]). Denote Zq = C the trivial representation 
and, for every j = 1, . . . ,n, define the following subspace: 

Zj = Y eL{Ln):f = f{xi,...,Xj), f{xi,...,Xj-i,i)=Oj 

of dimension r^i • • • mj-i{mj — 1). In virtue of the correspondence between 
^ui(r)-invariant operators and bi-5'to6^„i(7^)(0")-invariant functions, the 
corresponding eigenvalues are given by the spherical Fourier transform of 
the convolver that represents P, namely 

Xj = ^ P{xQ,x)4,j{x), 

where <j)j is the jth spherical function, for all j = 0, 1, . . . , 77. It is easy to 
verify that one gets: 

• Ao = 1; 

• Xj = 1 — aiQ2 • • • On-j , for every j = l,...,77— 1; 

• A„ = 0. 

In particular, if we set 

Pi = Qia2 • • • an-i(l - On-i+l) 

for every i = 1, . . . , 77 — 1, with p^ = l — ai and Pi = Ji for every z = 1, . . . , 77, 
the eigenspaces given for L(Xi x • • • x Xn) in Proposition 6.1 are exactly 
the Zj's with the corresponding eigenvalues. 

Let us prove that the eigenvalues that we have obtained in Proposi- 
tion 6.1 coincide with the eigenvalues corresponding to the eigenspaces 
Zq, Zi, . . . , Zn- 
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We want to get these eigenvalues by using the formulas given in Proposi- 
tion 6.1 for the eigenvalues of the nested product P by setting Pi = Ji, then 
p° = aia2 ■ ■ -an-iil — an-i+i) for i = 1, . . . , n — 1 and p^ = l — a\. First of 
all, we observe that the eigenvalues of the operator Pi are 1, with multi- 
plicity 1 and 0, with multiplicity mj — 1. So we get L{Xi) = Wq Wl, with 
dim(PF^^) = 77T,j — 1, for all i = 1, . . . ,n. Following the formulas that we have 
given, the eigenspaces of P are: 

. L{Xi) (g) LiX2) ® ■ • ■ L(X„_i) (8) Wl^; 

• L{Xi) (g) L{X2) (8) ■ • ■ L{Xn-j-i) W^-^ Wq~^^^ ®---®WS, for every 
j = l,...,n-l; 

The corresponding eigenvalues are: 

• Er=n-j+irf' foi' every j = 1, . . . , n - 1; 

We need to prove that, for every j = l,...,n — 1, the eigenvalue X^iLn-j+i 
is equal to 1 — 0102 • • • "j- We prove the assertion by induction on j. 

If J = 1, we have p^ = 1 — ai. Now suppose the assertion to be true for j 
and show that it holds also for j + 1- We get 

n n 

^ Pi = P° +Pn-j = 1 - 0102 • • -aj +ai • • -ajll - aj+i) 

i=n—j i=n—j+l 

= 1 — ai • • • ajOj+i. 

7. The second crested product. In this section we define a different kind 
of product of two spaces X and Y, that we will call the second crested 
product. In fact it contains, as particular cases, the crossed product and the 
nested product described in Sections 5 and 6, respectively. We will study a 
Markov chain P on the set 0^ of functions from X to y whose domains are 
fc-subsets of X, giving the spectrum and the relative eigenspaces. 

7.1. General theory. Let X be a finite set of cardinality n, say X = {1,2, 
. . . ,n}. For every k = 1, . . . ,n, denote by $7^ the set of /c-subsets of X, so 

that \nk\ = {i). 

Now let y be a finite set and let Q be a transition matrix on Y, which is in 
detailed balance with the strict probability r. Let Aq = 1, Ai, . . . , Xm be the 
distinct eigenvalues of Q and denote by Wj the corresponding eigenspaces, 
for every j = 0,1, . . . ,m, so that the following spectral decomposition holds: 

m 

L{Y) = ^W,. 

j=0 
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Moreover, assume that the dimension of the eigenspace associated with the 
eigenvalue 1 is 1 and set dim(VFj) = rrij, for every j = 1, . . . , m. 
Recall that the eigenspace Wq is generated by the vector 

(1^) 

\Y\ times 

and that Wj is orthogonal to Wq with respect to the scalar product {•,-)t, 
for every j = 1, . . . ,m. 

For every k = 1, . . . ,n, consider the space 

Gk = {{A,e):Ae Uk and 9 G Y^}, 

that is, the space of functions 9 whose domain A = dom{9) is a fc-subset of 
X and which take values in Y. 

The set = Ufc=o ©fc is a poset with respect to the relation C defined in 
the following way: 

C X if dom (if) C dom (x) and = x | dom{ip) ■ 

The Markov chain P on that we are going to define can be regarded as 
follows. Let <po <1 a real number. Then, starting from a function 9 E 0^, 
with probability po we can reach a function ip € 0^ having the same domain 
as 9 and that can differ from 9 at most in one image, according with the 
probability QonY. 

On the other hand, with probability 1 — po we can reach in one step a 
function G 0^. whose domain intersects the domain of in A; — 1 elements 
(on which the functions coincide), and in such a way that the image of the 
A;th element of the domain of (p is uniformly chosen. 

Note that P defines a Markovian operator on the space iv(0fc) of all 
complex functions defined on 0^. 

When Y is the ultrametric space, the Markov chain P represents the 
so-called multi-insect, which generalizes the Insect Markov chain already 
studied. In particular if \X\ =n, we consider k insects living in k different 
subtrees and moving only one per each step in such a way that their dis- 
tance is preserved, giving rise to a Markov chain on the space of all possible 
configurations of k insects having this property. 

In fact each element in 0^. can be regarded as a configuration of k insects 
and vice versa. For example, let 9 G 0^ be a function such that dom{9) = 
{xi, . . . , Xk} and 9{xi) = yi, with Xi G X and € y for alH = 1, . . . , k. Then 
the corresponding configuration of k insects has an insect at each leaf (xj, yi). 
They live in all different subtrees since Xi ^ xj for j. 

We observe that the cardinality of this space is (^)|y|*^. This space can 
be regarded as the variety of subtrees (see [4]) of branch indices {k, 1) in the 
rooted tree (n, 
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li 9,(p (z Qk, with domains A and B respectively, then define the matrix 
A, indexed by 0^, whose entries are 

1, ii \AnB\=k-l and 9\AnB = '^\AnB, 
0, otherwise. 

Observe that the matrix A is symmetric. 

The operator P can be expressed in terms of the operator associated with 
A and of another operator M as 

(6) P=poM + {l-po) 

norm[A) 

where M describes the situation in which the domain is not changed and 
only one of the images of the function 9 G 0^ is changed according with 
the probability Q on Y . An analytic expression for M will be presented 
below. On the other hand, A describes the situation in which we pass from 
a function whose domain is ^ to a function whose domain is ^ U {i} \ {j}, 
with i ^ A and j S A, and we choose uniformly the image in Y of the element 
i. So the action of A on 0,^ is an analogue of the Bernoulli-Laplace diffusion 
model. By norm{A) we indicate the number of nonzero entries in each row 
of the matrix associated with A. 

It is easy to check that M is in detailed balance with the strict probability 
measure defined as 

rM{e) = -^l[T{9{i)), 

where 9 G O^ and dom{9) = A. On the other hand, it follows from the def- 
inition of the Markov chain A that the weighted graph associated with A 
is connected. From this and from the fact that the nonzero entries of A are 
all equal to 1, we can deduce that A is reversible and in detailed balance 
with a uniform probability measure. This forces tm to be uniform and so 
we have to assume that r is uniform on Y and the matrix Q is symmetric. 

In this way, P is in detailed balance with the uniform probability measure 
vr such that Tr{9) = tttyItttx-) for every € Ofc. This choice of r guarantees that, 

if / is any function in Wj, with j = 1, . . . , m, then J2y£Y fiv) — 0- 

The spectral theory of the operator M has been studied in Section 5. 
In fact, it corresponds to choosing, with probability ^, only one element of 
the domain and to changing the corresponding image with respect to the 
probability Q on Y, fixing the remaining ones. So we focus our attention to 
investigate the spectral theory of the operator A. 
Let us introduce two differential operators. 
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Definition 7.1. (1) For every k = 2,.. .,n the operator : L(Gfc) — > 
L(0fc_i) is defined by 

for every F G L(6fc) and ip € Qk-i- 

(2) For k = 1, . . . ,n the operator — > L{Qk) is defined by 

for every F G L(6fc_i) and 6* G 9^. 

Observe that the operator D^. is the adjoint of Dj.. The following decom- 
position holds: 

l{@u) = l[y[ yA= L(y^). 

In order to get a basis for the space L{Y^), for every j4 G ri^, we introduce 
some special functions that we will call fundamental functions. 

Definition 7.2. Suppose that A^VLk and that F^ G L{Y) for every 
j G A. Suppose also that each F^ belongs to an eigenspace of Q and set 
fflj = \{j ^ A:F^ ^ Wi}\. Then the tensor product ^jgy^F-' will be called a 
fundamental function of type a = (ao, oi, . . . , am) in L{Y^). 

In other words, we have 

[^F\9) = \{P{e{j)), 

for every Q G Y^. We also set i{a) = ai + • • ■ + = ^ — oo- 

The introduction of the fundamental functions allows to give a useful 
expression for the operators M and A. 

If F G -L(y^) C L(efc) is the fundamental function F = (gj^g^F^', with 

\A\ = k and F^:Y — > C, then MF = i EjeA[((8)ieA,i^j ^0 ® Q-^^']- So, if 
9 G Gfc and dom{9) = A, we get 

Analogously one has {AF){9) = J2ipF{(p), where the sum is over all (p G 
0fc such that dom{(p) dom{9) = k — 1 and ip = 9 on dom{ip) n dom{9). If 
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A = {dom[9) n A) U {i} (we denote by U the disjoint union), then 

(a ((g) ) (^) = E P'i^) = n F'{m){Yl P'iy)] ■ 

Denote P^^a^A the subspace of L{y^^ spanned by the fundamental func- 
tions of type a and 

-Pfc,a = Pk,a,A- 

AeQk 



Lemma 7.3. Dk maps Pk^a to Pk-i,a', where g/ = (uq — l,ai, ... ,am)- 
Conversely maps Pk-i,a' to Pk,a- 

Proof. Let F be a fundamental function of type a in L(y^) and let 
B C A such that A = B U {i}. Then for every if € , we have 

e£Y^:<pcej<^A 

\yGY / jeB 

The value of J2yeY iu) is zero if G Wj for j = 1, . . . , m and so D^F = 
if ao = 0. If F' G VFo, then D^F € Pk-i,g/- 

Analogously, let F G Pk~i,a',B with i? G ^^fc-i- Then for every G Y^, 
A = B Li {i}, one has 

{DiFm= E 

= n ^'(^(j)) 

j<^B 

= F\e{i))l[F^{e{j)), 

where = 1 on y (and so G Wq). □ 

The restriction of to P^ a will be denoted by Dk „ and the restriction 
ofI)*toPfe_i,„, byl)^^,. 
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The study of the compositions of the operators Dk^a and Dl ^ plays a 
central role. In fact it will be shown that the eigenspaces of these operators 
are also eigenspaces for A. Consider, for example, -Dfc+i^^+i applied to a 
function F S L{Qk) and calculated on G 0^. The functions ip G Qk+i such 
that if^ 9 are in number of — k). Each of them covers k + 1 functions 

in Qk] one of them is the function 6, the other ones are functions in 0^ 
whose domains differ by the domain of 9 of an element and coincide on 
their intersection. These functions are in number of \Y\{n — k)k and they 
correspond to functions that one can reach starting from 9 in the Markov 
chain described by A. From this it follows that norm{A) = \Y\{n — k)k. 



Lemma 7.4. Let F G Pk,a,A, with A G 0^. Then 

Dl^Dk,a = \Y\{k-e{a))I + Qk,a, 
where Qk,a *s defined by setting 



(7) {Qk,aFm 



0, _ ifF'^Wo, 
\Y\Fi9), ifF'eWo, 



for every 9 £ Q^. such that \dom{9) D A\ = k — 1 and A \ dom{9) = {i}. We 
denote by 9 the function in 0^ whose domain is A and such that 9\A\{i} = 9 
and 9{i) = 0(io), where dom{9) \A = {«o}- 

Proof. Take F G Pk,a,A and 9 e@k- We have 

iDl,Dk,aFm= iDk,aF){ip) 
</jeefe_i:(/jce 

= E E ^H- 

dom{ijj)=A 

If dom[9) = A, then we get 

{Dl,Dk,aF){9) = Y.iT. iy)] n F\9{t)) 

= \Y\{k-l{a))'^F\9{t)) 

t^A 

= \Y\{k-e{a))F{9), 

where the second equality follows from the fact that J2y£YF-'{y) = \Y\ if 
F^ G Wo and J^yeY FHv) = whenever F^ ^ Wq. 
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On the other hand, if \dom{6) fl ^| = A; — 1, with A \ dom{9) = {i}, then 

{Dl^,D,^^F){e)=[Y,F\y)\ n F'm) 
Vyey / jeA\{i} 

0, _ iiF'i Wo, 
\Y\F{d), iiF'eWo, 

which is just the definition of Qk,a- D 

Lemma 7.5. Let F G Pk,a!,A, with A G $7^. Then 

Dk+l,aDl+l^a = \Y\{n - k)I + Qk,a, 

where Qk,a is defined as in (7). 

Proof. Take F € Pk,a',A and 6 e@k- We have 

{Dk+i,aDl^,^^F){e) = iD*k+i,aF){^) 

= E E ^M- 

dom{LL!)=A 

If dom{9) = A, then we get 

{Dk+i,aDl_,,^^F){9) = Y^i^) 

j&ACyeY 

= \Y\{n- k)F{9). 

On the other hand, if \dom{9) riA\=k — l, with A \ dom{9) = {i}, then 
{Dk+i,aDl^,,,F){9) = (y F\y)) J] ^'(^OO) 

0, _ ifF^^M^o, 
\Y\F{9), ifF^GTVo, 

This completes the proof. □ 

The fohowing corollary easily follows. 

Corollary 7.6. Let F e Pk,a',A, with A G (7^. Then 

Dk+i,aDl^,,a-Dla'Dk,a' = \Y\in + £(0) - 2k)L . 
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Consider now the operator Dk^a-Pk,a — ^ Pk-i,a'- 

Definition 7.7. For < £{a) <k<n, set 

Pk,a,k = Ker{Dk,a) 
and inductively, for k <h <n, set 

Ph,a,k = Dh,aPh-l,a',k- 

These spaces have a fundamental importance because they exactly con- 
stitute the eigenspaces of the operator A. This will be a consequence of the 
following proposition. 

Proposition 7.8. Ph,a'.k is an eigenspace for the operator Dh+i^aD'^^i ^ 
and the corresponding eigenvalue is \Y\{n + i{a) — k — h){h — k + 1) . 

Proof. We prove the assertion by induction on h.lih = k, from the last 
corollary we get Dk+l,aDl_^_^JP|^^,|^ = \Y\{n + e{a)-2k)I, since Dk,g/Pk,a',k = 
by definition of Pk,a',k- 

Now suppose the lemma to be true for k < t < h and recall that, by 
definition, we have Ph+i,a',k = Dh+i a'^h,g/' ,k- Moreover, Corollary 7.6 gives 

Dh+2,aDU2,a " ^/.+l,a' = \Y\{n + l{a) - 2{h + 1))/. 

So we get 

+ \Y\{n + £{a)-2ih + l))Ph+i^a',k 
= \Y\{n + i{a) -k- h){h -k + l)Dl^^^^,Ph,a",k 

+ \Y\{n + i{a)-2{h + l))Ph+i,a',k 
= \Y\{n + i{a) -k-h-l){h-k + 2)Ph+i,a',k, 

where the second equality follows from the inductive hypothesis and the 
third one from an easy computation. This completes the proof. □ 

Corollary 7.9. Ph. a' ^k o,n eigenspace for A of eigenvalue |l^|(n-l- 
£{a) -k- h){h -k + l) -\Y\{n - h). 

Proof. It suffices to observe that the operator Qh^a defined in (7) coin- 
cides with the operator A on the space Ph^a and then the assertion follows 
from Lemma 7.5 and Proposition 7.8. □ 
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In particular, after normalizing the matrix A we obtain norm{/\) ^^^^ 
the corresponding eigenvalue is |y|(„^_/t)/t (l^l(^ + ^{Ql} — k — h){h — k + \) — 
\Y\{n-h)). 

The following lemma holds. 

Lemma 7.10. Given i{a) and h, then, for i{a) <k< mm{h, "'+^(^) | ^ 
spaces Ph,a',k o-f^ mutually orthogonal. 

Proof. Each Ph, 

a' k is s-ii eigenspace for the self-adjoint operator Dh+i,a ^ 
^h+i a- Since the eigenvalue \Y\{n + i{a) — k — h){h — k + 1) is a strictly de- 
creasing function of k for k < ""'"^^-^ , then to different values of k correspond 
different eigenvalues. This proves the assertion. □ 

Recall that, if a = (oq, oi, . . . , am), we set a' = (oq — 1, ai, . . . , am) and, 
inductively, a^+^ = a'' - (1, 0, . . . , 0). 

Proposition 7.11. Let F be a function in P^a^^-k k- Then, for £(a) < 
k < '^"'"2^-^ and k <h <n + £(a) — k, we have 

,2 _ { n + £ia)-2k)\ih-k )\ 
{n + £{a) - k - h)\ 

In particular, a^h-i a' " ' ^k+i a''-''-i '^^ isomorphism of Pj^^^h^k f. onto 

Ph,a,k ■ 

Proof. We prove the assertion by induction on h. For h = k + 1 and 

F e Pk,a',k, we have 

= \Y\{n + iia)-2k)\\F\\^ 

by Proposition 7.8, so the assertion is true. For h> k+ 1, applying Propo- 
sition 7.8 to Di^ig.D'^a^ S6t 

= {Dh,aDh,aF>h~l,a' " " " -^fe+l.ah-fe-i -^^-l,a' " " ' F'l^i^g.'^-k-iF) 

= \Y\in + iia) -k-h + l){h- A:)||I)Li,a' • • • ■ 
Now the proposition follows by induction. □ 

Proposition 7.12. Assume i{a) <h< !i±M. Then: 
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(2) Dl_^-^^^ : Ph^a' — > Ph+i,a is an injective map. 

Proof. We prove the assertion by induction on h. 

Assume that (1) and (2) are true for i{a) -l<h<t< HlM::!, For 
h = i{a) — 1 we have Pe(^a)-i,a = so the proposition trivially holds. 

Since the operator D"^ ^ is the adjoint of D^ a we have the following de- 
composition: 

Ph,a = Ker{Dh,a) © I?;^,^^h-i,a' 

= Ph,a,h ffi „Ph-l,a' ■ 



In particular 
By induction 



Pt+l,a — Pt+l,a,t+l © P't+l,aPt,a'- 



t 

Pt,a' = ^ Pt,a',k 
k=i{a) 



and so 



Pt+l,a = Pt+l,a,t+l © Dt+l,a ( ® Pt,a',k 1 

\k=e{a) / 

t+1 

= ® Pt+l,a,k- 

k=e{a) 

This proves (1), while (2) follows from (1) and Proposition 7.11. □ 

Corollary 7.13. The dimension of the spaces Ph,a,k that appear in the 
decomposition of a is 



n + £{a) + l-2k fn\ f k \ f i{a 



n-k + 1 \k J \£{a) J \ai,...,a., 



Proof. From the previous proposition it follows that 
dim(Pt+i,a,t+i) = dim(Pt+i,a) - dim{Pt^a')- 

Now 

j = l 
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In fact, (j"]^) represents the number of {t + l)-subsets in X and (^^ J^^^ x 
njLi(dini(Vrj))°J represents the number of possible choices in the fundamen- 
tal function F = OrGA of functions belonging to the eigenspace Wi of 
L{Y). Thus 




Since, by Proposition 7.11, dim(Pft^a fc) = dim(P^ ^h-k f^) one can obtain the 
result replacing t hy k — 1 and a by a^~^. □ 

We want to find now the eigenvector of norm{A) associated with the eigen- 
value 1. Consider in -Pi,(i,o,...,o) the function 

n 

1=1 

where fi is the fundamental function of type (1,0,...,0) whose domain is 
{i}. Set 

(/) = -Pi,(i,o,...,o),o =: -^1,(1,0,. ..,o)-fo,(o,... ,0),0- 

So the element i^o = ^;j,(/,,o,...,o) ' " " ^3,{3,o,...,o)^2,{2,o,...,o)/ is the generator 
of the space Ph,{h,o,...,o),Oj which has dimension 1. Corollary 7.9 implies that 
Ph,{h,o,...,o),o is an eigenspace for ^^^^^^^ and the corresponding eigenvalue is 
1. Moreover, the connectedness of the graph associated with A implies that 
this is the unique (up to constant) eigenvector of eigenvalue 1. We denote 

by ^'i,(i,o,...,o),i the orthogonal subspace to Pi,(i,o,...,o),o in A,(i,o,...,o)- It has 
dimension n — 1. 

Observe that the definition of fundamental functions is strictly linked to 
the spectral theory of the operator Q and so of the operator M restricted 
to each domain. In fact, if F is a fundamental function in Ph,a.Aj with a = 
(ao, ai, . . . , Um) and A G Clh, then it is an eigenvector for the operator M and 
the corresponding eigenvalue is j^J2]^=o^j^j- So the set of the eigenvalues of 
M is given by (^) copies of these values. In particular, the eigenspace Ph,a,k 
of is also an eigenspace for M and an eigenvector in this space has 

eigenvalue j^J^J^o'^j^j- So, by Corollary 7.9 and definition (6) of P, we get 
the following theorem. 
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Theorem 7.14. Ph. a^k is an eigenspace for P with eigenvalue 

{a)-k-h){h-k + l)-{n-h) 



j=0 



h[n — h) 



Remark 7.15. It is easy to check that the operator M is not ergodic. 
In fact its associated graph contains (^) connected components and so the 
multiphcity of the eigenvalue 1 for M is (^) . 



On the other hand we aheady observed that the operator 



has 



norm(h.) 

the eigenvalue 1 with multiplicity 1. To conclude that it is ergodic it suffices 
to show that —1 is not an eigenvalue, that is, the associated graph is not 
bipartite. In fact consider 9 €Qh with domain {ii,...,ih} and 6{ij) =yj, 
for every j = 1, . . . ,h. By definition of A we can connect 9 with ip, whose 
domain is {ii, . . . ,ih-i,it}, ih 7^ H and such that ^{ij) = yj = 0{ij) for all 
j = 1, . . . ,h — I and f{it) = yt- Moreover 6 can also be connected with g 
whose domain is {ii, . . . ,ih-2,ihj'i't} and such that g{ij) = yj = 6{ij) for all 
j = 1, . . . ,h — 2,h and Q{it) = yt = ^{h)- On the other hand and g are 
connected as well and this proves that the graph is not bipartite. 

From Proposition 3.4 we can deduce the ergodicity for the operator P, 
since the multiplicity of the eigenvalue 1 is 1 and the eigenvalue — 1 does not 
appear in the spectrum of P. 

Remark 7.16. The second crested product reduces to the crossed prod- 
uct ii k = n and to the nested product if /c = 1. 

In fact, if k = n, the domain of a function 9 G 0„ cannot be changed 
and 6 can be identified with the n-tuple (yi, . . . , y^) G of its images. The 
operator P becomes 



1 



n ^ 



+1 



1 /« 



i=l 



which is the crossed product on the space y". 
If A; = 1, then A has the following expression: 



A: 



1 



Vi 



1 

0/ 



and norm{A) = n — 1. So we get 

P = Po{Ix(^Q) + {l-Po)(^ 



A 



norm{A) 



which is just the nested product of X and y , with Px 



^ and Py = Q. 
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7.2. Bi-insect. In what follows, we take y as a homogeneous rooted tree 
of degree q and depth m — 1 and we give an explicit description of the 
spectrum of the operator P = pqM + (1 — Po) norm{A) acting on the space 
L(G2). Therefore we are considering functions in 02 such that the image 
of each element of the domain is an insect. In other words, our aim is to 
diagonalize the bi-insect Markov chain defined in Section 2.2.2. 

Suppose X to be a set of cardinality n and let m > 3. Recall that we have 
the decomposition 

m—l 

L{Y) = W^, 



where Wq = C and 



/GL(L„_i):/ = /(xi,. 



q-l 
i=0 







for every j = 1, . . . , m — 1. Observe that dim(Wj) = q^~^{q — 1). 

The eigenspaces relative to the operator A/norm(A) are the subspaces 
of the form -P2,(ao,ai,...,am-i),fe' '^ith k = 0,1, 2. The corresponding eigenvalue 
is 

^ [g^'H^ + ^(«) - A; - 2) (2 - A; + 1) - q"'-\n - 2)]. 



^m— 1 J 



n 



2)2 

So, by dependence of 
• P 

^2,(ao,ai,...,am. 



(a), we get the following eigenspaces: 



0,2 



Po 



,ai, 



,<J7?l-l)ll 



ao = 0, with eigenvalue A : 

ao = 1, with eigenvalue A : 

ao = 2, with eigenvalue A : 

1, with eigenvalue A : 

2, with eigenvalue A : 



0, 



1 



2(n 
1 



n 



2' 



P 



2,{2,0,...,0),0 ^ '^0 



ao 
ao 
= 2 



n 



2{n-2) 



with eigenvalue A = 1. 



Now we describe the eigenvalues of these eigenspaces with respect to the 
operator M and to join the results. 

If F is a fundamental function of type (ao, ai, . . . , a^^i), then it has 
eigenvalue ^X^jl^^OjAj, where Xj = 1 — ^^r/_-^ is the eigenvalue of the 

eigenspace Wj, of dimension q-^~^{q — 1), occurring in the spectral decom- 
position of L(Y). From this we can fill the following table in which we give 
the eigenspaces, together with the corresponding eigenvalue and dimension. 
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(ao,ai,...,am-i),2- We have three different cases: 

(1) if oo = 0, the corresponding eigenspace is 

-f2,(0,...,0, ,o,...,o, ^^1^ ,0,...,0),2 

ith place jth place 

of dimension n(n — l){q — l)^q^~^q-^~^ , with eigenvalue ^(Aj + Aj); 

(2) if Oo = 1, the corresponding eigenspace is 

P2,{1,...,0, ,0,...,0),2 

ith place 

of dimension n{n — 2){q— l)q^~^, with eigenvalue Po^-^ + (1 ~Po) 2(n-2) ' 

(3) if uq = 2, the corresponding eigenspace is -P2,(2,o,...,o),2 of dimension 
with eigenvalue po + (1 — Po);^- 

,(ao,ai,...,am-i),i' havc two different cases: 

(1) if Oo = 1, the corresponding eigenspace is 

-f2,(l,...,0, ,o,...,o),i 

ith place 

of dimension n{q — with eigenvalue Pq^-i^ + 

(2) if ao = 2, the corresponding eigenspace is -P2,(2,o,...,o),i of dimension 
n - 1, with eigenvalue pQ + {l- Po) 2{n-2) ■ 

• -f2,(2,o,...,o),o- III this case, the dimension of the eigenspace is 1 with eigen- 
value 1. 
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